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Abstract 


The major sources of neutrons from plutonium dioxide nuclear fuel are con- 
sidered in detail. These sources include spontaneous fission of several of the Pu 
isotopes, (a, n) reactions with low Z impurities in the fuel, and (a, n ) reactions with 
18 0. For spontaneous fission neutrons a value of (1.95 ±0.07) X 10 3 n/s/g Pu0 2 
is obtained. 

The neutron yield from (a, n) reactions with oxygen is calculated by integrating 
the reaction rate equation over all a-particle energies and all center-of-mass angles. 
The results indicate a neutron emission rate of (1.14 ±0.26) X 10 4 n/s/g Pu0 2 . 

The neutron yield from («, n) reactions with low Z impurities in the fuel is pre- 
sented in tabular form for one part per million of each impurity. The total neutron 
yield due to the combined effects of all the impurities depends upon the fractional 
weight concentration of each impurity. The total neutron flux emitted from a par- 
ticular fuel geometry is estimated by adding the neutron yield due to the induced 
fission to the other neutron sources. 
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Neutron Radiation Characteristics of 
Plutonium Dioxide Fuel 


I. Introduction 

One method of providing electrical power to spacecraft 
for outer planet missions is to utilize the heat produced 
in the radioactive a-decay of 238 Pu in a Pu0 2 -fueled 
radioisotope thermoelectric generator (RTG). Significant 
amounts of radiation are also given off by the fuel, includ- 
ing neutrons, y-rays, and gamma radiation accompanying 
the a-decays. This paper is concerned with accurately 
characterizing the neutron radiation in particular. 

In order to obtain the neutron yield from Pu0 2 fuel, it 
is necessary to determine the number of neutrons resulting 
from each of the two major sources, spontaneous fission of 
the plutonium isotopes, which produce a “natural” neutron 
background, and (a, n ) reactions with light elements 
present in the fuel. These light elements include oxygen as 
well as trace quantities of numerous other impurities con- 
tained in actual fuels. All of these neutrons are capable 
of causing induced fissioning in the plutonium, giving rise 
to still another source of neutrons. A fourth source of neu- 
trons is from (y, n) reactions, although, as will be shown, 
this source is ' not significant. The neutron-producing 
reactions can therefore be summarized as (1) spontaneous 
fission, (2) induced fission, (3) (a, n) reactions with low Z 
isotopes, and (4) photoneutron formation. 

In order to obtain the spectral distribution of the neu- 
tron radiation from a Pu0 2 heat source, the effects of all 


the above reactions must be combined. In this report each 
of these neutron spectra will be investigated individually 
and then combined to yield the composite spectrum. To 
simplify the calculation of these neutron spectra, any of a 
number of schemes may be adopted to eliminate one or 
more of the sources. For example, because of large 
surface-to-volume ratio, induced fission effects may be 
eliminated by choosing a disk-like source of Pu0 2 . More- 
over, by utilizing a plutonium metal source, such as the 
HP-15-2 source used in various phases of artificial heart 
research, only spontaneous fission neutrons will be of 
concern. 


II. Spontaneous Neutron Radiation from 
238 Pu Isotopes 

The maximum binding energy per nucleon occurs at 
an atomic number of about 60. In heavier nuclei the total 
binding energy can be increased by dividing the original 
nucleus into two smaller nuclei. These spontaneous fission 
reactions do not occur in the common elements because 
they are opposed by a potential barrier. In heavy ele- 
ments such as 238 Pu, the dissociation energy is above the 
barrier height, and spontaneous fission, in which neutrons 
are emitted, can take place. A plot of partial half-life for 
sponstaneous fission against the fission parameter Z 2 /A 
indicates that odd-A nuclei have a much longer half-life 
for spontaneous fission than do even-A nuclei with the 
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same fission parameter. In the case of 238 Pu and 239 Pu, the 
half-lives for spontaneous fission are about 5 X 10 10 and 
10 16 years, respectively. 1 

These fissionable elements, however, are also a-emitters 
with much shorter half-lives (i.e., 86.4 years for 238 Pu and 
2.4 X 10 4 years for 239 Pu). Thus one expects the spontane- 
ous radiation to be less important than (a, n) reactions. 
Nevertheless, one cannot ignore the fact that these spon- 
taneous fission neutrons are capable of induced fission and 
can have energies as much as 10 MeV or higher. This 
means they can cause considerable damage to scientific 
packages if they are exposed for long periods of time. 


Integration of this equation over all energies gives the 
total neutron yield from 238 Pu spontaneous fission. A frac- 
tion of the neutrons with energy greater than 5meV is 
due to the induced fission neutrons, but the shape of the 
spectrum of these neutrons is the same as for the spon- 
taneous neutrons (Fig. 1). 


Using Eq. (1), the average neutron yield may be evalu- 
ated; that is, 


y = 



dE = 2.8 X 10 3 


(2) 


in n/s/g 238 Pu. 


A. Spontaneous Neutron Spectra of 238 Pu 

Since plutonium metal used as a heat source and fur- 
nished by AEC does not emit pure spontaneous fission 
neutrons, the mixed radiation of the plutonium heat source 
must be analyzed. Clearly, it is possible to obtain the neu- 
tron fission spectrum of a plutonium power source by 
fitting a Maxwellian distribution to the experimental data 
and thereby obtaining the required parameters to simu- 
late the neutron spectrum (Ref. 1). However, since a pure 
plutonium source is difficult to make because of the lack 
of techniques for removal of impurities, the contribution 
from a mixture of impurities and different plutonium iso- 
topes will cause the fitted parameter to be incorrect or 
biased. Thus a scheme which allows the least bias of im- 
purities and mixture of isotopes has to be used. 

Anderson and Neff (Ref. 2) of Monsanto Research Cor- 
poration used a HP15-2 metal radiation power source to 
obtain the pure spontaneous fission neutrons. The neutron 
spectrum of this source was composed of neutrons from 
both spontaneous fission of 238 Pu and from (a, n) reactions 
with impurities in the fuel. The researchers estimated that 
75% of these neutrons were from spontaneous fission and 
that these made up essentially all of the spectrum for 
energies greater than 5 MeV. Thus by using this argument 
and data points from 5 to 14 MeV, they used a weighted 
least-square fit to a Maxwellian distribution in order to 
approximate the spontaneous fission spectrum. The result 
of this fit is found to be: 

N(E) = 2.04 X 10 3 (E)% exp ( - E/1.34) (1) 

in n/s/MeV/g 238 Pu, where E is the energy of the emitted 
neutrons in MeV. 


'A more accurate value of spontaneous fission half life of “Pu has 
been obtained by Los Alamos and Argonne Laboratories. 


Since the maximum energy of neutrons from 18 0 (a, n) 
is about 4.5 meV and only 12% of the neutrons were esti- 
mated to be due to (a, n) reactions with impurities, 2 the 


2 From a nonrelativistic two-body reaction and by using 5.5-MeV« 
particle energy (Ref. 3). 



Fig. 1. Neutron spontaneous fission spectrum 
(n/s/MeV/g 238 Pu) 
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assumption that beyond 5 MeV nearly all the neutrons are 
from spontaneous emission was justified and, therefore 
Eq. (1) may be used as the spontaneous fission spectrum 
equation. 

A Maxwellian distribution is assumed here because the 
neutrons emitted from the nucleus are modeled to a Fermi 
gas of temperature T, and therefore the normalized energy 
distribution is given by r) N ( E ) : 

Vir ( E ) = 2 ir* T-% EK exp ( — E/T) (3) 

where T is the corresponding nuclear gas temperature. 

This equation depends only on the one parameter T 
which, on comparison with Eq. (1), would be equal to 
1.34 MeV for spontaneous emission of neutrons from 238 Pu. 
The upper limit of the errors on this parameter from all 
sources had been estimated (Ref. 2) to be about 5%. Utiliz- 
ing this error, the theory and experiment set a value of 
(1.34 ±0.07) Mev for the nuclear temperature of 238 Pu. 

The average energy is calculated by using the nor- 
malized energy spectrum and is equal to (2.0 ±0.1) MeV; 
the most probable neutron energy corresponds to the peak 
value of the neutron energy distribution, E p = T/2 = 
(0.67 ±0.03) MeV. 

B. Spontaneous Fission Parameters of Other 

Fuel Components 

So far, the spontaneous fission parameters of neutrons 
from 23S Pu have been calculated. By having the Max- 
wellian distribution one is able to obtain neutron yield 
and average energy, or by using the average number of 
neutrons per fission one is able to calculate the energy 
spectrum, if the fission parameter is known. 


Since a pure 238 Pu isotope in a fuel furnished by AEC 
as a heat source does not exist, the contribution from other 
plutonium isotopes mixed up with the 238 Pu should be 
considered. If a similar nuclear emission model is assumed 
for each of the isotopes, one needs only two parameters 
to estimate their emission rates or energy spectra: the 
fission parameter v and the neutron yield Y. Table 1 sum- 
marizes these parameters for the spontaneous fission of Pu 
isotopes; Table 2 presents the spontaneous fission yield of 
all the plutonium isotopes in the fuel. 

C. Combined Spontaneous Neutron Source Spectrum 

In order to obtain the combined spontaneous neutron 
emission from a nuclear fuel consisting of a number of 
isotopes, the radiation spectra of all components should 
be combined according to their fractional weight and half- 
life ratio. Weightwise, the only other component isotope 
which may contribute a significant amount of radiation is 
23B Pu, but the ratio of half-life of 239 Pu/ 238 Pu is equal to 
about 10 -5 . Therefore, it is clear that only 238 Pu need be 
used to obtain spontaneous neutron spectrum and yield. 
This is emphasized by the fact that the total contribution 
of neutrons from spontaneous emission is only a fraction 
of the total neutron yield. 

In Table 1, various values of fission parameters such as 
number of neutrons/fission v, average energy of emerging 
neutrons E, the most probable energy E p , half-life 
and the corresponding values of nuclear temperature are 
quoted for different plutonium isotopes included in the 
power source. Even though these values are within 1 035 
of their average values, the radiation spectrum of each 
isotope is calculated separately because of tremendous 
differences in their half-lives and percent weight contri- 
bution to the overall weight of the fuel. For 238 Pu isotope, 
v = (2.33 ±0.08) neutron /fission, f — 1.165 X 10 3 fission/ 
g/s, the neutron yield y = (2.71 ±0.09) X 10 3 n/g 238 Pu/s, 


Table 1. Spontaneous fission parameters for plutonium isotopes 
(Furnished by AEC as a heat source) 


Pu isotopes 

Abundance", % 

v (n/tl" 

E, MeV 

E p , MeV 

T, MeV 

Power, W/g 

236 

1.2 X 10' 4 

2.30 ±0.19 

1.96 

mam 


18 

238 

80.27 

2.33 ±0.08 

1.965 



0.55 

239 

15.87 

2.88 ±0.016 

2.065 

0.69 

1.38 

1.9 X 10' 3 

240 

3.022 

2.257 ±0.046 

1.955 

0.655 

1.31 

6.9 X 10 J 

242 

0.643 

3.03 ±0.03 

2.08 

0.695 

1.38 

4.7 X 10' 3 

242 

0.132 

2.18 ±0.09 

1.946 


1.30 

1.2 X lO’ 4 


“Obtained from Ref. 4. 
b Obtained from Ref. 5. 
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Table 2. Plutonium spontaneous neutron fission yield 

Pu isotopes (i) 

Abundance, % 

T yr 

(1) 

fi, 

f T /T* 

*238* V4 ' * % 

( 2 ) 

n/g PuVs 

Neutron yield 
(3) 

n/g PuVs 

(41 

n/g PuOj/s 

236 

1.2 X 10 -* 

3.5 X 10 ” 

1.43 X 10 ” 

3.3 X 10 4 

3.93 X 10~ 2 

3.46 X 10- 2 

. 238 

80.27 

4.3 X 10 “ 

1.16 X 10 ” 

2.71 X 10 3 

2.18 X 10 s 

1.92 X 10 * 

239 

15.87 

5.5 X 10 “ 

9.1 X 10 -’ 

2.63 X 10~ 2 

4.19 X 10 J 

3.60 X lO "* 

240 

3.02 

1.22 X 10 " 

4.1 X 10 ” 

9.27 X 10 2 

2.80 X 10 1 

2.46 X 10 1 

241 

0.643 

5.0 X 10 “ 

1.0 x io- J 

3.04 X 10- 2 

1.95 X 10- 4 

1.77 X 10- 4 



(estimated) 





242 

0.143 

6.8 X 10 “ 

7.87 X 10 ” 

1.61 X 10 s 

2.3 X 10 ° 

2.02 X 10 ” 





Total 

2.21 X 10 ” 

1.95 X 10 * 


(1) : the number of fissions/g PuVs. 

(2) : the neutron yield for Pu*. 

(3) : the neutron yield for Pu* in the fuel. Table 1. 

(4) : the neutron yield for each plutonium isotope dioxide in the fuel. 


and for a 2200-W (thermal) metal fuel, disk-shaped power 
source, M = 3.52 X 10 3 g, the total neutron rate would be 
(9.60 ±0.34); X 10 6 neutron/s. Multiplying the neutron 
yield by the fraction of 238 Pu in Pu, namely 0.81, a total 
of 2.2 X 10 3 n/g Pu/s is obtained. Moreover, if the neu- 
tron yield for the spontaneous fission per gram Pu per sec- 
ond is multiplied by the fraction of Pu in Pu0 2 , (0.8814), 
the neutron yield from spontaneous fission will be equal 
to 1.95 X 10 3 n/g Pu0 2 /s. 

III. Analytical Calculation of the Neutron Yield 
from ( a , n) Reactions with Low-Z Impurities 

Impurities are usually found in a plutonium heat source 
because of the lack of techniques to remove them com- 
pletely. Clearly, (o, n) reactions can take place with all 
impurity elements in the fuel when the energy of the 
incoming a particle is greater than the threshold energy 
of the reaction. However, both the height of the potential 
barrier of the target nucleus and the a-ray energy should 
be within an allowable range. Therefore, it may be as- 
sumed that the neutron yield from high-Z impurities is 
small, and the element with atomic number equal to 14 
(i.e., silicon) is the highest element for which the neutron 
yield need be calculated (Ref. 6). 

A. Impurity Characteristics 

The a particle emitted from the plutonium isotope has 
various energies with different probabilities of emission. 
The energy and the emission rate of the a particle from 
the 238 Pu decay process depend upon the excited states 


of the daughter nucleus, namely, 234 U isotope. The 238 Pu 
isotope decays into the ground state of 234 U isotope 71% 
of the time; the corresponding a particle has an energy 
of 5.491 MeV (Refs. 7-9). If the daughter nucleus is left 
in the first excited state, the a-particle energy will be 
5,448 MeV, with an emission rate of 29%. Table 3 sum- 
marizes these values for the 238 Pu to 234 U decay scheme. 
The emission rates from other excited states of 234 U are 
negligible, and thus an average energy of 5.48 MeV can 
be used for the a-particle energy emitted from 238 Pu. An 
a particle, when passing through an impurity element 
or plutonium isotopes, loses energy by ionization and 
excitation processes. This means the energy loss of a 
particles plays an important part in the neutron yield 
determination. 

The y-ray energy emitted from the excited states of 
the daughter nucleus ( 234 U*) has a maximum energy of 


Table 3. Alpha-emission characteristics of 
238 Pu -» 234 U* -» 234 U decay scheme 

ay 

(Ref. 8) 


Excited 
state of 

Ground 

state 

First 

excited 

state 

Second 

excited 

state 

Third 

excited 

state 

Fourth 

excited 

state 

Fifth 

excited 

state 

a 

emission, 

% 

71 

29 

0.13 

0.005 

7 X 10 -” 

10- 4 

MeV 

5.491 

5.448 

5.35 

5.2 

5.0 

— 

£,* 

MeV 

0 

0.0435 

0.293 

0.43 

0.499 

0.81 
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43.5 keV and thus will be absorbed by the source or the 
shield surrounding it. Photons are also radiated from dif- 
ferent plutonium isotopes in the fuel and from the (a, n ) 
reactions; these are not the subject of this study. 

Table 4 summarizes different parameters in the a-decay 
process (Ref. 10) of all the important plutonium isotopes 
considered. The energy of the a particle was averaged 
over the probabilities of a emission, and Fa is calculated 
to indicate the effect of other isotopes besides 238 Pu on 
the neutron spectrum in general. The F„ has been evalu- 
ated by multiplying the ratio of half-lives by the percent 
weight of each plutonium isotope. From the values of the 


last column, it can be concluded that the contributions to 
the total number of («, n) reactions are mainly due to the 
238 Pu isotope. 

Table 5 summarizes the list of low-Z impurities which 
have an acceptable threshold energy for (a, n) reactions. 
The potential barrier height B is equal to 

Zz e 2 (r 0 A% + 1.2)- 1 

where e 1 — 1.44MeV-Fermi, r 0 = 1.3 Fermi, Z and z are 
the target and a-particle atomic number, respectively, and 
r 0 is the nuclear radius. In the classical sense, the (a, n) 


Table 4. Characteristics of Pu isotopes as a-emitters 


Pu isotopes 

Weight percent 

a-activity, a/ g/s 

f % <«>« Y' 

a, P, or y 
emission 

E weighted, 
MeV 

ratio 

F„ 

236 

1.2 X 10 * 

1.96 X 10" 

2.85 

a 

5.73 

33 

4.0 X 10-“ 

238 

80.27 

6.4 X 10” 

86.4 

(89.6) 

a.y 

5.51 

1.0 

8.027 X 10’ 1 

239 

15.87 

2.28 X 10’ 

2.44 X 10* 

a.y 

5.51 

3.6 X 10” 

5.70 X 10 " 

240 

3.022 

8.43 X 10* 

6.58 X 10* 

a 

5.18 

1.32 X 10-’ 

3.96 X 10” 

241 

0.643 

1.7 X 10’ 

13.2 

0.0024 a 
fS (99%) 

4.893 

6.54 

4.20 X 10” 

242 

0.132 

1.47 X 10’ 

3.79 X 10’ 

a 

4.87 

2.79 X 10-’ 

3.68 X 10” 


Table 5. (a, n) reaction with impurities” 

Impurity 

z 

A 

(m — A) 

X 10 5 , amu 

(E.)th, 

MeV 

(£n)innx, 

MeV 

Q, 

MeV 

B, 

MeV 

Abundance, % 
(Ref. 1 1 ) 

Reaction 

li 

3 

7 

18.223 

4.39 

3.57 

-2.79 

2.33 

92.1 

Li 7 (a, n) B 10 


3 

8 

25.018 

— 

13.54 

6.63 

2.27 

~0.0 

Li’(a,n)B" 

Be 

4 

9 

15.043 

— 

12.47 

5.71 

2.95 

100.0 

Be’ (a, n) C” 


4 

10 

16.711 

— 

10.45 

3.85 

2.88 

~0.0 

Be 70 (a, n)C" 

B 

5 

10 

16.114 

— 

7.55 

1.07 

3.60 

18.4 

B'”(<», n) N” 


5 

11 

12.789 

— 

6.51 

0.15 

3.52 

81.6 

B” (a, n) N” 


5 

12 

18.162 

— 

14.17 

7.36 

3.45 

0.0 

B” (a, n) N“ 

C 

6 

13 

7.473 

— 

8.54 

2.20 

4.06 

1.1 

C" (a, n) O” 


6 

14 

7.515 

2.55 

4.07 

— 1.98 

3.99 

0.0 

C" (a, n) O” 

O 

8 

17 

4.533 

— 

6.71 

0.60 

5.07 

0.04 

O” (a, n) Ne” 


8 

18 

4.857 

0.86 

5.31 

-0.70 

5.00 

0.204 

0“ (a, n) Ne" 

F 

9 

18 

6.651 

3.12 

3.34 

— 2.56 

5.63 

0.0 

F” (a, n) No" 


9 

19 

4.456 

2.33 

4.00 

— 1.92 

5.55 

100.0 

F” (a, n) Na" 


9 

20 

6.350 

— 

10.14 

3.90 

5.48 

0.0 

F 70 (a, n ) Na" 

Na 

11 

23 

— 2.945 

— 

9.16 

3.15 

6.40 

100.0 

Na" (a, n) Al" 

Mg 

12 

25 

—6.255 

— 

8.52 

2.67 

6.90 

11.5 

Mg" (a, n) Si" 


12 

26 

9.190 

— 

5.93 

3.99 

6.84 

11.1 

Mg" [a, n) Si" 


12 

27 

-7.124 

— 

10.19 

4.22 

6.78 

0.0 

Mg" ( a,n ) Si” 

Al 

13 

28 

—9.230 

— 

7.88 

1.96 

7.27 

100.0 

Al" (a, n) P" 

Si 

14 

29 

-14.340 

1.80 

4.22 

-1.58 

7.76 

6.2 

Si" (a, n) S” 


14 

30 

-16.763 

4.22 

1.94 

-3.73 

7.70 

4.2 

Si” (a, n) S” 


14 

31 

— 14.850 

— 

7.05 

1.18 

7.63 

0.0 

Si” (a, n) S" 

“Atomic mass numbers are obtained from Ref. 12. 
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reaction is not allowed if E a < B even though E a < E a ; 
however, from quantum mechanics only the value of the 
transmission coefficient can determine the transparency of 
the Coulomb potential barrier (CPB). The transmission 
coefficient for the S wave a particle is T 0 = e~ r , where y 
is the “Gamow” factor and is equal to 



where /? = v„/C. Since (E a ) max = 5.5 MeV, (/? mln = 0, 
/?ma* = 0.05) and (y min = +1.83Z, y max = go). This means 
the transmission coefficient has a range of 

T 0 = [0, exp ( — 1.83Z)] 

For silicon, T 0 is 10 -11 while for Cd it is 10 -39 ; therefore, 
one can neglect the contribution from those impurities 
which have a very small probability of penetrating the 
potential barrier. 


B. (a n) Neutron Yield 

The neutron yield from impurities y i in a mixture can 
be written in terms of the experimentally measured neu- 
tron yield from target impurity atoms, the activity of the 
plutonium radioactive nuclei, and the straggling ratio; 
that is (Refs. 3, 6, 13), 


Yi K (\N) Pu y 0 . ^ + SpuNpu ^ (4) 

where K is the fraction of the a-emitter, (kN) Pu is called 
the a-activity of the emitter source, and Y 0j is the experi- 
mentally measured neutron yield for the impurity in ques- 
tion. The last fraction in the expression above accounts 
for the effect of competition between the (a, n) reactions 
and other unwanted ionization processes by which a par- 
ticles also lose energy and is termed the straggling ratio. 
Therefore, the neutron source density depends upon the 
degree of mixing of the ^-particle source and impurities. 
The term X is the disintegration constant and is equal to 
0.693 /T w ; N is the total number of radioactive atoms of 
the emitter. The term (XZV) is given in units of curies (one 
curie is equal to 3.7 X 10 10 disintegrations per second) 
and X.N = —dN ( t)/dt , where N ( t)/N (0) = exp ( — A#) is 
the number of radioactive atoms left at time t. The term 
S 4 is the ith element atomic stopping power and is propor- 
tional (Ref. 6) to Z{ (Z { + 7)~ % . 


Table 6 presents the experimental values of the n/a 
yield ( y 0 ) for different low-Z materials, as reported by 
different groups. The results of Gorshkov (Ref. 14), owing 
to its detailed experimental study and small errors re- 
ported, are employed here. Equation (4) is used to evalu- 
ate neutron yields for all the impurities involved, and the 
yields are compared with other measured data. 


Table 6. Observed neutron yield for low-Z impurities' 1 
(Ea = 5.298 Mev from Polonium) 


Element 

z 


Neutron yield, n/10 

6 a 

A 

Gorshkov (Ref. 14) 

Roberts 
(Ref. 11) 
(corrected) 

Breen, Hertz, 
and Wright 
(Ref. 15) 

Li 

3 

7 

— 

2.6 

1.0 

Be 

4 

9 

84.4 ±0.9 

80 

59 

B 

5 

10 

19.6 ±0.2 

24 

15 

C 

6 

12 

0.113 ±0.015 

0.11 


o 

8 

16 

0.068 ±0.011 

0.07 

— 

F 

9 

19 

11.6 ±0.2 

12 

— 

No 

11 

23 

— 

1.5 

1.1 

Mg 

12 

25 

1.33 ±0.04 

1.4 

0.80 


13 

27 

0.16 ±0.03 

0.14 

0.53 

Si 

14 

28 

16.5 ±0.02 

16 

— 

Cl 

17 

35 

— 

0.11 

— 

A 

18 

36 

— 

0.38 

— 

P 

15 

31 

— 

— 

— 

S 

16 

32 

— 

— 

— 

"A is the atomic mass number for the most abundant isotope in an 
element (Table 5), and yields are for a mixed isotope target. 


In order to simplify Eq. (4), it may be assumed (Ref. 13) 
that (SN)i < < (SN) Pu , and since Avogadro’s number N 0 
is equal to ( NA/p ) (Ref. 3) 


t/i = K\ Pu y 0 . 


Si TTijiVo 

Spu A» 


yi 

m Pu 


Si 


KXpuf/oj g 


m j N n 
m Pu Ai 


(5) 



( 6 ) 


where Yj is the neutron yield of the ith impurity per gram 
of plutonium per part impurity, and i/ 0j is the observed 
atomic neutron yield of the ith impurity and is given 
in neutron yield /a-particle. Thus the total neutron yield 
of an impurity is given by Eq. (6) in neutrons/gram of 
plutonium/second/part of impurity. 
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C. Total Neutron Yield 

The observed number of neutrons per a-particle in the 
(a, n) reaction with impurities has been reported by differ- 
ent groups. Table 6 shows the observed atomic neutron 
yield in n/10 6 <* for impurities up to silicon. 

Through the use of Eq. (6) and Tables 4 and 6, neutron 
yields of all impurities capable of producing neutrons can 
be obtained. By knowing the neutron yield of each im- 
purity, the total neutron yield of the actual fuel (where 
the a-emitter is mixed with other impurities) can be cal- 
culated. This is achieved simply by adding the contribu- 
tions from all the impurities to obtain the total yield, 
namely, 

Y = )(£),, (7) 


isotope is used. The half-life of 238 Pu is 86.4 yr, and 1 yr is 
3.15 X 10 7 s. Therefore, 

0,693 

^ Ty, (second) 

= 2.54 X 10 -10 (second) -1 

\N 0 = 1.52 X 10 14 second -1 mole -1 

If one takes K as 81* 238 Pu in the fuel, Eq. (8) can be 
rewritten in the form of Y (N/g 238 Pu/s/ppm): 

i 

for 1 ppm of the ith impurity; i.e., fi = l. Table 7 sum- 
marizes the computer results obtained from this equation. 


where / { is the fractional weight of the ith impurity, In the calculation of neutron yield, the experimental 
K\ Pu depends on the a-emitter, y 0 is the measured num- values obtained by Gorshkov (Ref. 14) were used. The 

ber of neutrons/a, and S f /S Pu depends on the Z-number absolute errors on the neutron yields given in Table 7 are 

of each impurity. Since S, is proportional to Z (Z + 7) - ^ about 1.5*. This is due to the small errors given by 

Gorshkov. 

The data in Table 7 indicate that the light elements 
(8) such as Be, B, and F must be eliminated or minimized in 
the fuel to preclude large neutron emission rates. It should 
be noted that although the yield from oxygen per ppm is 
Through the use of Eq. (6) and Tables 4 and 6, the neu- relatively low, the percent concentration of this element 

tron yield for one gram of plutonium 238 and one part in Pu0 2 fuel is substantial and will tend to contribute sig- 

per million (1 ppm) of each impurity has been calculated nificantly to the overall neutron yield. Clearly, for a given 

with the aid of a computer code; the results are pre- fuel source the neutron yield of each impurity should be 

sented in Table 7. Since only the a particles emitted by calculated according to its contribution in parts per mil- 

238 Pu are of concern, the disintegration constant A. of this lion to the overall fuel concentration. In the last column 

Table 7. Neutron yield due to (a, n) reactions with impurities in a 23S Pu heat source 11 


[Matlack and Metz (Ref. 6)], 




Element 

z 

A 

n/10 8 a 

Si 

(?). 

r, 

n/s/g 333 Pu/ppm 

Sp« 

This work 

Matlack and 
Metz (Ref. 6) 

Li 

3 

7 

2.6 ± 1 

0.10 

3.72 X 10 -1 

4.64 ±1.78 

4.6 

Be 

4 

9 

84.4 ±0.9 

0.13 

9.35 

148.80 ±1.55 

133 

B 

5 

11 

19.6 ±0.2 

0.15 

1.78 

32.80 ±0.33 

41 

C 

6 

13 

0.113 ±0.015 

0.18 

8.7 X 10 -s 

0.20 ±0.03 

0.2 

O 

8 

16 

0.065 ±0.011 

0.20 

4.27 X 10' 3 

0.10 ±0.01 

0.1 

F 

9 

19 

11.6 ±0.2 

0.24 

6.10 X 10"' 

18.07 ±0.31 

18 

Na 

11 

23 

1.5 ±0.4 

0.28 

6.55 X 10" 3 

2.22 ±0.58 

2.2 

Mg 

12 

25 

1.33 ±0.04 

0.29 

5.33 X 10 -3 

1.93 ±0.06 

2.1 

At 

13 

28 

0.76 ±0.03 

0.31 

2.82 X 10 -3 

1.08 ±0.04 

1.0 

Si 

14 

29 

0.168 ±0.02 

0.32 

5.8 X 10 -3 

0.24 ±0.03 

0.2 

P 

— 

— 

— 

— 

— 

— 

< 0.03 

S 

— 

— 

— 

— 

— 

— 

< 0.03 

“The definition of parameter A and the value of B a are given in Table 6. 
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of Table 7, the results reported by Matlack and Metz 
(Ref. 6) are presented. Other experimental work in this 
area has been performed by E. D. Arnold (Ref. 16) and 
T. W. Bonner, A. A. Kraus, Jr., J. B. Morion, and J. P. 
Schiffer (Ref. 17). 


IV, Neutron Yield from (a, n) Reactions with 
ls O Isotope 

When the fuel is made up of Pu0 2 , the oxygen isotope 
plays an important part in the emission of neutrons via 
the ls O (a, n) 21 Ne reaction. The analysis of this interaction 
is given elsewhere (Refs. 13, 18, 19). In this section, the 
result of such a calculation is analyzed in terms of the 
total neutron yield. 

In a hypothetical experiment it was assumed that the 
alpha particles from the 238 Pu isotope are allowed to inter- 
act with the oxygen isotope in the 238 Pu0 2 compound 
used as a nuclear power source. When an alpha particle 
enters into a nucleus, its energy dissipates among all 
nucleons. The excess energy brought in by the alpha par- 
ticle will cause one or more nucleons to come to the sur- 
face and escape. The excited nucleus then emits a photon 
and returns to the ground state. 

A 5.48-MeV alpha particle is used at the start of the 
(a, n) reaction, and it is assumed that the energy loss takes 
place by ionization interaction with the target atoms. A 
computer code was developed (Ref. 19) to integrate 
Eq. (10) and to obtain the total neutron yield at a given 
neutron energy E„(k ): 


Yi [E n ( k )] = 2.66 X 10- 2 — XJa 



Ea<j mb ( E a ) (cos B* k y — cos 6* k) dEa 
8 In (5.62 Ea) + 47 In (0.547 E a ) 


(10) 


where N is the number density, Y a is the a-yield of the 
238 Pu isotope, and X t is the probability of the decay of the 
ith excited state of the 21 Ne recoil nucleus (Ref. 20). 


Since N T /N 0X = 1.815 X 10 3 , Eq. (10) may be rewritten 
in the form of 


Yi [E„ ( k )] = 4.82 X 10- 10 X<Y J, ( E„, max ; E*, J (11) 

where /< is the value of the integral given in Eq. (10). 
Table 8 presents different values of the decay probability, 
alpha particle threshold energy, Q-value, and, finally, the 
maximum energy of the emitted neutrons for each excited 
state of the recoil nucleus. In Eq. (11), the neutron yield 
Yj and the a-yield Y a are given in N/s/g 238 Pu and a/ s/g 


Table 8. Excited states of the recoil isotope 
in ls O (a, n) 21 Ne reaction 


Excited states of 
**Ne isotope 

E, MeV 

Probability 
(Ref. 20), 
% 

Q, MeV 

(Ea)th, 

MeV 

MeV 

Ground state 

0.0 

55 

—0.704 

0.86 

4.48 

First excited 

0.349 

35 

— 1.049 

1.28 

4.11 

state 

Second excited 

1.75 

~10 

— 2.45 

3.00 

2.60 

state 

Third excited 

2.87 

~1 

—3.57 

4.36 

1.33 


state 

238 Pu, respectively. Therefore, the number of neutrons per 
alpha per MeV from the ith excited state of the 21 Ne iso- 
tope produced in 18 0 (n, a) 2l Ne reaction is 

,i = 4.82X 10- 10 XiI i (E (l ) (12) 

Figure 2 presents different values of the function rji in 
terms of the emitted neutron energy E„ for each excited 



0 1 2 3 4 5 6 7 

NEUTRON ENERGY, MeV 


Fig. 2. Weighted neutron yield from 18 Ola, n) 21 Ne 
reactions 
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state of the recoil. The area under each curve represents 
the number of neutrons emitted per a particle interacted. 
The sum of the areas gives the total neutron yield for the 
reaction. 3 From this analysis one may obtain the neutron 
yihld for the (a, n) reaction in a Pu0 2 nuclear power source 
(Fig. 2 and Ref. 18). 

As far as the experimental values are concerned, con- 
sider first the shape of the yield spectrum. One such 
experiment was performed at the Savannah River Labora- 
tory (SRL, Ref. 21). Figure 3 presents the result of the JPL 
calculation and adjusted measured data (shown as the 
dotted curve). One can easily notice a discrepancy be- 
tween the two in the low-energy range, namely, below a 
1-MeV neutron energy. It has been explained (Ref. 21) 
that these low-energy neutrons are due to the scattering 
and multiplication in the fuel itself. 

In order to investigate the source of these neutrons the 
SRL neutron source geometry and the FASTER Monte 
Carlo code (Ref. 22) have been used to obtain the fraction 
of neutrons scattered in the source. Figure 4 shows three 
neutron fluxes, all calculated by using the Monte Carlo 
code and the computed neutron source. A significant dif- 
ference can be noted when the unscattered flux is com- 
pared with the scattered flux obtained when neutrons are 
allowed to scatter five times. It can be seen from Fig. 4 

3 For more detailed discussions, Refs. 18 and 19 should be consulted. 



Fig. 3. Neutron yield from the SRL experiment 


that if the number is allowed to increase to 15 scatterings 
no major difference is observed. By use of the unscattered 
flux, the fraction of neutrons scattered as a function of 
energy is obtained. This scattering percent can now be 
folded into the computed neutron yield and compared 
with the experimental value shown in Fig. 3. In the com- 
parison, note that the high-intensity, low-energy neutrons 
cannot be accounted for by the scattering of the neutrons 
in the fuel. 

In addition, the average path length of neutrons can 
be estimated to be about 2 cm, and the average energy 
change due to each neutron scattering is about 10%. Since 
the diameter of the fuel geometry used by SRL is also 
2 cm, one cannot expect a large number of neutrons to be 
scattered and lose a lot of energy; this probably cannot 
account for the high intensity of neutrons in the low- 
energy range. 

One may account for these low-energy neutrons as 
being due to the thermal neutrons which are, for example. 



Fig. 4. Monte Carlo result calculated for the SRL 
geometry with JPL neutron yield spectrum 
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created around the experimental apparatus. This is made 
more pronounced by the fact that the detector selected for 
this experiment was a lithium detector which has a high 
(3-value; thus it is overbiased toward low-energy neutrons. 

In order to compare the integrated value of the neutron 
yield spectrum with the experimental data, one has to 
obtain the area under each curve (Fig. 2) and then sum 
them all up. 

Figure 2 presents the weighted neutron yield ( n/a ) from 
each of the excited states of the 21 Ne-isotope considered. 
The total integrated values for each of these states are 
also shown in Fig. 2; they are 9.55 n/10 6 <* for the 21 Ne 
ground state and 6.07, 1.62, and 0.094 n/10 6 a for the first 
three excited states. Thus a total of (17.33 ±4.0) n/10 6 a 
is obtained from 238 Pu 18 0 2 when the contributions from all 
the states are added together. 

The magnitude of the error on the total neutron yield 
is due to the errors on the reaction cross section mea- 
sured by Bair and Willard (Ref. 23). Their value for the 
cross section at the peak of the 3.19 MeV resonance is 
(109 ±25) mb. 

It can be noted that usually the measurements have 
been performed on the Pu0 2 fuel, where Pu contains all 
of the plutonium isotopes of which 81% is the 238 Pu iso- 
tope. Since the fractional weight of Pu in Pu0 2 is equal 
to 0.8814, the fractional weight of 238 Pu in Pu0 2 is equal 
to 0.714. The fraction of ls O in the oxygen is 1.815 X 10~ 3 ; 
thus the neutron yield from the Pu0 2 source is equal to 
(2.24 ±0.51) X 10- 8 n/a. 


Table 9. Neutron yield comparison 
in the ls O (a, n) 21 Ne reaction 


Source 

n/s/g “Pu, 
10‘ 

n/s/g Pu, 
!0 4 

n/s/g PUO 2 , 
10 4 

This work 




Y„ - 6.35 X 10" 
(Ref. 6) 

This work 

(1.99 ±0.45) 

(1.67 ±0.38) 

(1.42 ±0.32) 

y„ = 5.1 1 x io" 
(Ref. 9) 

(1.60 ±0.37) 

(1.34 ±0.31) 

(1.14 ±0.26) 

Matlock and Metz 
(Ref. 6) 

. 

— 

(1.00 ±0.02) 

Anderson and Neff 
(Ref. 2) 

(2-3) 

(2.2 ±0.1) 

— 

Stoddard and 

Albensius 
(Ref. 9) 

— 

2.1 

1.9 


10 


The number of a particles per gram of plutonium 238 
per second, the o-yield, is (Ref. 6) 6.35 X 10 11 ; therefore, 
the total neutron yield from the (a, n) interaction with 
oxygen is (1.42 ±0.32) X 10“ n/g Pu0 2 /s. An experimental 
value given by the Los Alamos group is (1.0 ±0.02) X 10 4 
n/g Pu0 2 /s which is in good agreement with the com- 
puted result given here. 

Moreover, if a value of 5.11 X 10 11 is used (Ref. 9) for 
the a-yield of 238 Pu, the calculated value will be equal to 
(1.14 ±0.26) X 10 4 n/gPu0 2 /s, and the agreement will be 
even better. These results are presented in Table 9. 

V. The Overall Neutron Yield and Flux 
from PuO, Fuel 

In the previous sections some of the schemes which 
may be used to obtain the partial neutron yield and yield 
spectrum from a Pu0 2 fuel power source were indicated. 
The physical phenomena behind the emission of the neu- 
trons from such a source are the spontaneous fission, in- 
duced fission, and (a, n ) reactions with impurities in the 
fuel. 

The only natural neutron radiation of this fuel source 
is the spontaneous fission; however, the presence of im- 
purities in the fuel makes the problem more complex. 
Presented in this section is a way to add the neutron yield 
spectra of all these sources together and, for a typical 
multihundred-watt power source, determine the neutron 
flux at various observation points outside the geometry 
of the source. 

A. Total Neutron Source 

The total neutron source for a given fuel depends upon 
the percent impurities in the fuel. If the fuel is free of any 
impurities which contribute to the total neutron yield and 
if the 18 0- 17 0 isotopes are eliminated or reduced to a 
negligible amount, the total neutron yield spectrum is 
identical to the spontaneous fission yield. If all the con- 
tributing impurities are removed and the ls O isotope is 
left in the fuel, the total neutron yield spectrum will be 
equal to the sum of the spontaneous fission yield and the 
(a, n) reaction with oxygen. Such a result is presented in 
Table 10, where the computed results and the measured 
values reported by the Los Alamos group (Ref. 24) are also 
shown. 

In the case where not all other additive impurities are 
removed, one has to know the exact ppm contamination of 
these impurities. Since only the low-Z impurities need be 
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Table 10. Neutron yield from a Pu0 2 fuel power source 
(n/g Pu0 2 /s) 


Source* 1 

Spontaneous 

fission, 

10* 

la, n) with 
oxygen, 
10 4 

Total, 

10 4 

JPL: Ya f® Pu) = 
6.35 x ro 11 
(Ref. 6) 

(1.96 ±0.10) 

(1.42 ±0.32) 

(1.62 ±0.37) 

JPL ; y„ f®Pu) = 
5.1 1 X 10" 
(Ref. 9) 

(1.96 ±0.10) 

(1.14 ±0.26) 

(1.34 ±0.31) 

Los Alamos 
experimental 
(Ref. 24) 

(2.20 ±0.04) 

(1.00 ±0.02) 

(1.22 ±0.02) 

“The T„ is fhe 

Pu a yield given in a / g 

"'Pu/s. 


considered (Section III), the problem becomes somewhat 
simplified; that is, the high-Z impurities do not contribute 
to the overall neutron yield from a Pu0 2 power source 


unless they contain oxygen. In such cases the oxygen por- 
tion of the impurity can be treated as an impurity, and the 
neutron yield contribution can be evaluated via (<*, n) 
reactions with oxygen. 

In addition, the neutron yield from impurities emitted 
by spontaneous fission or («, n) reactions, where the a par- 
ticle is emitted from the impurities, should not be con- 
sidered as a major source, basically because of a very low 
percent contamination of these undesirable impurities in 
the fuel. The enhancement of the total neutron yield from 
these sources is much less than 1%. In order to calculate 
the total neutron yield, including the contribution from 
the impurities, a typical multihundred-watt (MHW) fuel 
capsule with the contaminating impurities shown in 
Table 11 has been selected. Using the neutron yield/ppm 


Table 1 1 . Fractions of low-Z impurities in a 
2.2-kW(thermal) Pu0 2 fuel power source 


Impurity 

elements 

Impurities in 
PuOz, ppm 

Yield/ppm, 
n/g PuOj/s/ppm 

Yield 

Li 

i 

4.64 

4.64 

Be 

10 

148.80 

1488.00 

B 

20 

32.80 

656.00 

C 

500 

0.20 

100.00 

O 

11750* 

0.10 

1175.00 

F 

250 

18.70 

4517.50 

Na 

300 

2.00 

600.00 

Mg 

50 

1.93 

96.50 

Al 

400 

1.08 

432.00 

Si 

441 

0.240 

105.84 

Others 

1200 

0.03 

36.00 




9211.50 ±138.2 


‘Oxygen content of the fuel other than the PuOs. 


for each impurity, one can estimate the total integrated 
neutron yield for this source at slightly less than 10 4 n/g 
Pu0 2 /s. If the neutron yield spectrum due to the ( a,n ) 
reaction with oxygen is used for the neutron yield spec- 
trum due to the (a, n) reactions with other impurities, the 
total neutron yield spectrum can be computed by adding 
the three spectra together. The result of such a computa- 
tion is shown in Fig. 5. 

The neutron yields per gram of Pu0 2 per second for 
spontaneous fission and (a, n) reaction with oxygen and 
other impurities are shown in Fig. 5. The total neutron 
yield spectrum is then obtained by adding the neutron 
spectra mentioned above. The average neutron energy is 
equal to 2.48 MeV. When the induced fission neutron and 
scattering within the neutron source are included, this 
average energy is reduced to 2.28 MeV. For any other fuel 
power source design, the ppm impurities have to be 
evaluated first; a new neutron yield due to the impurities 
can then be easily calculated from Table 11. 



Fig. 5. Total neutron yield from a 
2.2-kW-MHW P u0 2 source 
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B. Neutron Flux from a MHW Fuel Capsule 

The total neutron yield spectrum obtained in the pre- 
vious sections for a given source can now be used to 
calculate the neutron flux. As was pointed out previously, 
a 2200-W (thermal) MHW fuel capsule was chosen for 
this calculation as a typical neutron source. The geometry 
of such a source is shown in Fig. 6. The dimensions of the 
source and shielding materials are also shown. The fuel 
volume is about 707 cm 3 , the power density is 3.11 W 
(thermal)/cm 3 , and the specific power is 0.39 W (ther- 
mal)/g Pu0 2 . The two areas of observation are at 50 and 
100 cm away from the center of the source, and in each 
segment three points are used. The point at the center is 
45 deg apart from the other two points at the opposite 
comers. 

Points 1 and 2 are located right at the outside of the 
shielding material and about 2 cm away from the source 
itself. These two points are the nearest points to the source 
for which the flux calculation is made. 



Fig. 6. The geometry setup of the 2.2-kW-MHW Pu0 2 
fuel capsule 


1. Neutron multiplication. The real neutron source is 
different from an extremely thin disk-like hypothetical 
source because induced fission and scattering can occur. 
These induced fission neutrons are created in the same 
fashion that the spontaneous fission neutrons are emitted. 
Therefore, additional neutrons are emitted via induced 
fission, and the overall neutron yield is increased by a 
multiplication factor (MF). Because the shape of the in- 
duced fission spectrum is the same as the spontaneous 
fission spectrum, one may simply add them together. 

If Y t is the total integrated neutron yield when the 
induced fission is included and Y 0 when the induced fission 
Y if is not included, then Y t is given by 

Y ( -Y 0 + Y„ (13) 

and 

Y 0 = Y(«,n) + Y, 

For the ith neutron energy interval the corresponding 
relationship will be 

Y t (i) = [(Y* Y- 1 - 1) Y„ Y-] Y f (i) + Y„ (i) (14) 

where Y/ and Y it are the neutron yields for the sponta- 
neous and the induced fission, respectively. Using Eq. (14) 
one can obtain the total neutron yield spectrum including 
the induced fission neutrons. The only parameter which 
is not known is Y,Y“ 1 or the multiplication factor which 
has to be calculated by utilization of a transport code. 
Such a computer code has been developed by Atomics 
International and is called ANISN (Ref. 25). It is used to 
obtain a diffusion solution of the transport equation 
utilizing a first-order iteration technique. 

After feeding the geometry of the MHW neutron source 
and Y 0 , the neutron yield, into the ANISN code, a multi- 
plication factor of 1.81 was calculated for this source. 
The reason for such a high multiplication factor has to do 
with the number density of the fuel, the surface-to-volume 
ratio, and the quantity of graphite used for shielding 
purposes. Graphite causes the slowing down of neutrons 
and consequently increases the induced fission neutron 
yield. For this configuration of the fuel geometry the 
corresponding K e f/ is equal to 0.448. This means the in- 
duced fission is about 44.8% of the total neutron yield. If 
the graphite is removed, the ANISN code gives a K et f of 
0.353, which corresponds to a multiplication factor of 1.55. 
This factor is about 20% larger than the VIKING SNAP 
19 multiplication factor obtained from the ANISN code. 
The SNAP 19, of course, has only 675 W(thermal), in com- 
parison with the 2200 W and a much smaller density and 
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Table 12. Total neutron yield for a 2.2-kW(thermal) MHW PuO., fuel power source 


Neutron 

Neutron yield, 
n/s/g Pu0 2 


Neutron yield spectrum, 
n/s/g PuOa/MeV 

MeV 

Total 

Total 

effective® 

Total 

Total 

effective® 

(2.05 — 10) X 10~ 8 

2.53 X 10- 8 

2.90 X 10"’ 

3.36 X 10-’ 

3.37 X 10° 

(1.00 — 30) X 10-’ 

4.21 X 10" 6 

4.85 X 10-” 

1.45 X 10° 

1.67 X 10’ 

(3.00 — 30) X 10 ° 

1.37 X 10- 4 

1.58 X 10“ 3 

5.07 X 10° 

5.84 X 10’ 

(3.00 — 55.5) X 10 5 

2.12 X 10- 2 

1.34 X 10”’ 

4.04 X 10’ 

2.55 X 10 2 

(5.55 — 170) X 10- 4 

2.83 X 10° 

2.14 X 10’ 

1.72 X 10 2 

1.30 X 10 3 

(1.70 — 4.49) X 10 ” 

1.14 X 10’ 

8.08 X 10’ 

4.08 X 10 s 

2.90 X 10 3 

0.0449-0.122 

6.22 X 10’ 

3.93 X 10 s 

8.07 X !0 2 

5.09 X 10 3 

0.122-0.201 

8.29 X 10’ 

4.97 X 10 2 

1.05 X 10 3 

6.30 X 10 3 

0.201-0.331 

1.70 X 10 2 

9.93 X 10 2 

1.31 X 10 3 

7.64 X 10 3 

0.331-0.546 

3.51 X 10 2 

1.85 X 10’ 

1.63 X 10 3 

8.58 X 10 3 

0.546-0.702 

3.08 X 10 2 

1.46 X 10 3 

1.98 X 10 3 

9.33 X 10 3 

0.702-0.900 

5.09 X 10 2 

1.95 X 10 3 

2.57 X 10 3 

9.35 X 10 3 

0.900-1.16 

8.84 X 10 2 

2.70 X 10 3 

3.40 X 10 3 

1.04 X 10 4 

1.16-1.49 

1.52 X 10 3 

3.60 X 10 3 

4.59 X 10 3 

1.09 X 10 4 

1.49-1.91 

2.69 X 10 3 

4.98 X 10 3 

6.41 X 10 1 

1.19 X 10* 

1.91-2.45 

5.28 X 10 s 

7.63 X 10 3 

9.78 X 10 3 

1.41 X 10 4 

2.45-3.14 

7.67 X 10 3 

9.84 X 10 3 

1.11 X 10* 

1.43 X 10 4 

3.14-4.04 

5.13 X 10 s 

6.94 X 10 3 

5.70 X 10 3 

7.71 X 10 3 

4.04-4.46 

5.26 X 10 2 

1.07 X 10 3 

1.25 X 10 3 

2.55 X 10 3 

4.46-5.18 

6.30 X 10' 

7.26 X 10 s 

8.75 X 10’ 

1.01 X 10 3 

5.18-6.66 

6.86 X 10’ 

7.90 X 10 2 

4.64 X 10’ 

5.34 X 10 2 

6.66-8.55 

3.23 X 10’ 

3.72 X 10 2 

1.71 X 10’ 

1.97 X 10 2 

8.55-10.00 

6.95 X 10° 

8.01 X 10’ 

4.79 X 10° 

5.52 X 10’ 

Total 

2.54 X 10 4 

4.59 X !0 4 



Average £„ 

2.48 

2.28 



“When the induced fission is included. 


volume/surface ratio. It should be emphasized that the 
ANISN code is one-dimensional, and consequently the 
margin of error is substantial. 

2. Flux calculation. By using the procedure described 
above, the total neutron yield and yield spectrum were 
calculated for the 2200- W MHW fuel capsule and are 
presented in Table 12. Having obtained the total neutron 
yield from the MHW fuel power source, the Monte Carlo 
code now can be applied to calculate the neutron flux at 
vari bus observation points shown in Fig. 6. The FASTER 
Monte Carlo (Ref. 22) code was selected for this com- 
putation. 

In Table 13 the total number and energy flux of the 
selected detector points are shown. In addition, in Fig. 7 
the flux spectrum of the detector is presented. Clearly, 
for points 1 and 2 (i.e., closest points to the shield) the 
1/R 2 falloff hypothesis does not hold; for the points much 
farther away from the source, the 1/R 2 falloff may be 
utilized to obtain the flux. Thus the integrated values of 
neutron flux need to be calculated only for the smaller 
region surrounding the source (i.e., about 50 cm). For the 


Table 13. Neutron flux at nine detector points 


Detector 

points 

<*r 

cm 

8, 

deg 

<f>n 

n/cmVs 

<f>E, 

MeV/cm 2 /s 

Bn, 

MeV 

i 

6 

0 

1.25 X 10“ 

2.56 X 10 5 

2.06 

2 

10 

90 

4.83 X 10 4 

8.75 X 10 4 

2.06 

3 

50 

0 

5.43 X 10 3 

1.08 X 10‘ 

2.00 

4 

50 

45 

3.73 X 10 3 

7.55 X 10 3 

2.02 

5 

50 

90 

1.45 X 10" 

3.10 X 10 3 

2.14 

6 

100 

0 

1.39 X 10 s 

2.77 X 10 3 

2.00 

7 

100 

45 

9.10 X 10 2 

1.84 X 10 2 

2.03 

8 

100 

90 

2.62 X 10 2 

5.70-X 10 2 _. 

2.18 - 

9 

25 

45 

1.80 X 10 4 

3.60 X 10* 

2.00 


approximate values of the neutron flux beyond these 
points (first five points), 1/R 2 falloff can be easily applied 
(Ref. 19). 


C. Results and Conclusion 

The total integrated neutron yield and spectrum have 
been calculated for a MHW fuel power source with in- 
duced fission included. The multiplication factor was 
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E„, MeV 

Fig. 7. Neutron flux spectrum at 50-cm away from 
the 2.2-kW Pu0 2 fuel capsule 

calculated by using the ANISN transport code. The total 
neutron yield spectrum and the physical geometry of a 
2.2-kW(thermal) MHW source were fed into the Monte 
Carlo code to obtain the neutron flux and the neutron flux 
spectrum at various points outside the fuel. It is concluded 


that, beyond 50 cm, a point source assumption can be 
used. For neutron fluxes at these large distances the 1/R 2 
falloff can be used with confidence. 

Calculations have indicated a total of 2.54 X 10 4 n/s/g 
Pu0 2 when the induced fission neutrons are not included. 
This yield consists of 1.95 X 10 3 n/s/g Pu0 2 due to the 
spontaneous fission, 1.42 X 10 4 n/s/g Pu0 2 due to the 
(a, n) reactions with oxygen in Pu0 2 , and 9.21 X 10 3 n/s/g 
Pu0 2 due to the ( a , n) reactions with all other impurities. 
The induced fission neutron yield for this typical fuel 
power source is equal to 2.05 X 10 4 n/s/g Pu0 2 . There- 
fore/by adding the induced fission neutron yield to the 
total neutron yield, a grand total of 4.59 X 10 4 n/s/g Pu0 2 
is obtained. In order to obtain these neutron yields in 
terms of 238 Pu rather than Pu0 2 , the above numbers 
should be multiplied by 1/0.714. 

From these data it is concluded that 4.2% of the neutrons 
are emitted because of spontaneous fission, 31% via the 
(a, n) reaction with oxygen in Pu0 2 , 20% via (a, n ) reactions 
with impurities in this fuel configuration, and, finally, 
44.8% via induced fission. The K e ,, is calculated to be equal 
to 0.448, which gives a multiplication factor of 1.81. 

The total neutron rate from a 2.2-kW-MHW source 
may be obtained by multiplying the neutron yield 
(4.59 X 10 4 n/s/g PuO z ) by the total mass of the Pu0 2 fuel 
(2200/0.39 = 5.65 kg). The neutron rate is thus equal to 
2.60 X 10 s n/s. Moreover, for 2-MeV neutrons a 1 mrem/h 
dose rate is produced by 8.7 n/cm 2 /s (dose rate conver- 
sion factor); therefore, at any detector point (Table 12), 
the total dose rate may be calculated. For example, at 
100 cm distance from the source (i.e., detector point 7) 
the number flux is 910 n/cm 2 /s, and the total dose rate at 
this location is 104 mrem/h. 
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